Quantum Error Correcting Codes From The Compression Formalism 
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Introduction. — Quantum computers will rely on a smor- 
gasbord of error correction techniques to combat harmful 
effects such as decoherence on bits of information that 
are physically encoded within quantum systems. The 
quantum error correction "toolbox" now includes sev- 
eral strategies for accomplishing these feats, but there are 
still many deep issues to resolve. The standard method 
for error correction via active intervention in the quan- 
tum computing regime has been cleanly phrased in terms 
of an analysis of operators on Hilbert space. For every 
quantum channel that arises through the usual system- 
environment formalism, there is a family of error (or 
"Kraus" ) operators that describe the possible corruption 
by the channel of qubits encoded as states in, or oper- 
ators on, the system Hilbert space. The main protocol 
for quantum error correction (QEC) E 0. 13 depends 
upon the existence and identification of states and opera- 
tors on which the error operators are jointly well-behaved 
in a precise sense. 

The stabilizer formalism for QEC 0, gives a con- 
structive framework to find correctable codes for error 
models of "Pauli type". While there are other success- 
ful techniques that can be applied in special cases (for 
instance, see 0, II H E3, El El El El), the landscape 
of general strategies to find codes for other classes of 
channels is fairly sparse. In particular, the theory lacks 
a systematic method that applies to arbitrary quantum 
channels. Indeed, after spending any time at all on this 
problem, it becomes clear that this is an extremely daunt- 
ing challenge. 

Nevertheless, in this paper, we introduce an approach 
based on a "compression formalism" that may lead to 
such a general method. In particular, we cast the general 
problem of finding correctable codes for quantum chan- 
nels into a matrix analysis framework. We then utilize a 
new tool recently introduced in E3 ~ called the "higher- 
rank numerical range" - the study of which was primarily 
motivated by this problem. As an application, we solve 



the quantum error correction problem in its entirety for 
the class of "bi-unitary channels" on two-qubit (i.e., four- 
dimensional) Hilbert space, and construct qubit codes for 
such channels in the arbitrary dimension case. In the case 
of Pauli-error models, we show that this approach cap- 
tures codes not obtained via the Pauli matrix stabilizer 
formalism. 

The rest of the paper is organized as follows. Next we 
recall basic properties of quantum channels and the for- 
mulation of the active quantum error correction protocol. 
We follow this by introducing the higher-rank numerical 
range as a tool in quantum error correction. We then 
consider randomized unitary channels, and focus on the 
bi-unitary case. This is followed by a characterization 
of correctable codes in the two-qubit case, and a consid- 
eration of the connection with the stabilizer formalism. 
We discuss possible further applications and limitations 
of this approach throughout the paper, and finish with a 
conclusion. 

Quantum Channels and Active Error Correction. 
Consider an open quantum system S represented on a 
Hilbert space Ti, and write B(Ti) for the set of operators 
that act on Ti. A "snapshot" of a Hamiltonian- induced 
evolution of S is called a quantum channel. Mathe- 
matically, channels are represented by completely pos- 
itive, trace preserving maps £ on B{TL). (For experimen- 
tal reasons, the current focus in quantum computing is 
on finite-dimensional Hilbert spaces, and thus we shall 
make this assumption throughout the paper.) The struc- 
ture theorem [16L Il7| for completely positive maps shows 
that every quantum channel £ on Ti has an operator- 
sum representation of the form £(o~) = ^ Q E a oE\ for 
all a S B(Ti), where the "error" operators (or "Kraus" 
operators) E a are operators that act on Ti. A set of 
error operators {E a } can always be chosen with cardi- 
nality at most N 2 := dim(7i) 2 . For simplicity, we shall 
always assume the maximum cardinality holds, by possi- 
bly including zero operators as some of the E a . The trace 
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preservation condition is equivalent to ^2 a E\E a — J. As 
a notational convenience, we shall write £ = {E a } when 
the E a determine £ through the operator-sum represen- 
tation, and, as a further convenience, wc will also use 
this notation when scalar multiples of the E a are error 
operators for £. 

In the context of quantum computing, the operators 
E a are the errors induced by the channel £. Thus, error 
correction protocols in quantum computing are crafted 
primarily to mitigate the effects of such operators on 
quantum information encoded in evolving systems. By 
"active quantum error correction" , we mean protocols 
that involve active intervention into the system to cor- 
rect errors. The basic method for active quantum error 
correction 0, 0, El Q identifies quantum codes with sub- 
spaces C of the system Hilbert space TL. Then a code C 
is correctable for a channel £ if all states encoded in C 
prior to the action of £ can be fully recovered in a man- 
ner allowed by quantum mechanics. From the operator 
perspective, C is correctable for £ if there is a quantum 
channel 1Z on TL such that 



(lZo£)(a) = a, 



(1) 



for all operators a supported on C; that is, all a of the 
form a = PqoPq where Pq is the projection of TL onto C. 

There is a very useful characterization 0, 3 of cor- 
rectable codes for a given quantum channel £, when a 
set of error operators £ = {E a } is known. The code C is 
correctable for £ = {E a } if and only if there is a scalar 
matrix A = (A a &) such that 



Pr.EtE b P c = \ nh Pr. Va, b 



(2) 



Thus, Eq. (J2J shows how the physical problem of active 
quantum error correction can be phrased in terms of a 
clean mathematical statement that involves relations sat- 
isfied by the error operators. The prototypical scenario 
occurs when C is correctable and the matrix A = Ac is 
diagonal. Here the error operators take the code space 
to mutually orthogonal subspaces E a C. In the case that 
each of the errors restricted to C is a scalar multiple 
of a unitary U a , the correction operation 1Z is given by 
1Z = {UlP a }, where P a is the projection of TL onto the 
subspace U a C. 

Let us note that every matrix Ac which satisfies Eq. 
for some code space C is necessarily positive. Indeed, 
define E to be the N 2 x N 2 positive block matrix 



E 



[Eh. E 2 ■ • •] > 0, 



(3) 



where Ei,E 2 , ... is (any) enumeration of the set {E a }, 
and [Ei E 2 ■ • ■ ] is the operator row matrix mapping from 
H (n2) to TL. (So EjEj is the entry of E.) Then 



observe that the set of equations from Eq. J5J can be 
succinctly stated as the single matrix equation 



< P C EP C = A c ® Pe, 



(4) 



where Pc is the N 2 x N 2 diagonal block matrix with Pq 
as the operator in each of the diagonal entries. In fact 
Ac is a density matrix when the E a satisfy the trace 
preservation constraint (and C is non-zero); 

(E A -) P C = J2 P C E aEaP C (5) 

a a 

= P c (J2ElE a )P c = P c . (6) 



Higher-Rank Numerical Range and Projections. — The 
characterization from Eq. J2J of correctable codes moti- 
vates consideration of the following notion. Let a belong 
to B{TL). For k > 1, define the rank-k numerical range 
of a to be the subset of the complex plane given by 

Afc(cr) = {A e C : Po-P = XP for some P S V k ), (7) 

where Vk is the set of all rank-fc projections on TC. We 
refer to elements of Afe(cr) as "compression- values" for er, 
since they are obtained through compressions of a to k- 
dimensional subspaces. The case k = 1 yields the familiar 
numerical range W(o~) for operators [18j; 



Ai(a)=W{o-) = {(aiP\iP) 
It is clear that 



G H, 



A 1 (a)DA 2 (a)D...DA N (a). 



1}. (8) 



(9) 



Of course, the cases k > 1 are of immediate interest 
in quantum error correction. We refer to the sets Afc(cr), 
k > 1, as the higher-rank numerical ranges. This no- 
tion was recently introduced in where a number of 
mathematical properties were developed. Here we briefly 
outline the main points. The following facts from |l5| ap- 
ply to arbitrary a £ B(H). 

Proposition 1 Let a E B(C N ). For all k > 1, the rank- 
k numerical range Afc(cr) is a compact subset of the com- 
plex plane C. If 2k > N, then Afe(cr) is either empty, or 
a singleton set. If Afc(er) = {Ao} is a singleton set with 
2k > N, then Ao is an eigenvalue of geometric multiplic- 
ity at least 2k — N. In particular, An (a) is non-empty if 
and only if a is a scalar matrix. 

The main result of 0] is stated as follows. 

Theorem 1 If a = <r' is a Hermitian operator on 7i — 
C N , with eigenvalues (counting multiplicities) given by 



at < a 2 < • • ■ < aw, 



(10) 



and k > 1 is a fixed positive integer, then the rank-k nu- 
merical range Afe(cr) coincides with the set [a k , ajsr-k+i] 
which is: 
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(i) a non-degenerate closed interval if ak < ajy^k+i , 

(ii) a singleton set if = O/v-fe+ij 

(iii) an empty set if > ajy-k+i ■ 

The cases of non-degenerate spectra for TV = 4 and 
N = 6 are depicted in Figure 1. We note that the proof 
from 0| of this result is constructive, in the sense that 
the projections P associated with each A e Afc(cr) can be 
explicitly identified. We shall extend this fact for special 
cases in the analysis of the next sections. 
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FIG. 1: Numerical range Afc(cr) for a non-degenerate Hermi- 
tian operator a of size a) N = 4 and b) N = 6 with spectrum 
{a*}. 



Let us consider the following simple examples for the 
purpose of illumination. 

(1) Let o~ = o~ be the operator on 2-qubit space H = 
C 2 ® C 2 given by 

a = |00)(00| + 2|01)(01| + 3|10)(10| + 4|11)(11|. 

Then Ax(cr) = [1,4], A 2 (cr) = [2,3], and A 3 (cr) and 
A^(a) are empty. 

(2) Consider the Pauli matrix Z = |0)(0| - |1)(1|. Then 
Ki{Z) = [-1,1] and A 2 (Z) is empty. On W = 
C 2 <g> C 2 , if we let Z x = Z H 2 , then Ai(Zi) = 
[-1,1], A 2 (Zi) = [-1,1], and A 3 (Zi) and A 4 (Zi) 
are empty. More generally, if Z\ = Z % l 2 «-i acts 
onM= (C 2 )®", then A fc (Zi) = [-1,1] for k < 
2 n ~ 1 . 

We have extended TheorernQJin 0| partly to the case 
of normal operators; i.e., a S £>(7i) such that aa^ = cr^er. 
Recall that the convex hull coT of a subset T C C is 
given by 

m 

coT= |*iAi + ...+t m A m : y^Jj = l,tj > 0, 

3=1 

Aj e r,m > i|. 



Notice that the following result applies to unitary op- 
erators. We have conjectured in [15j that the converse 
inclusion of Eq. 1)11(1 below holds for arbitrary normal 
operators a. 

Theorem 2 If a is a normal operator onH = C , then 
for all k > 1, i/ie rank-k numerical range Afe(cr) is a 
subset of every convex hull co F where T is an (JV+1 — fe)- 
point subset (counting multiplicities) of the spectrum of 
a; that is, 



AfcO) c p| ( co r). 



(11) 



In the next section we shall present a solution to the 
error correction problem for a special class of channels. 
The solution depends on the analysis of higher-rank nu- 
merical ranges discussed above and specific properties of 
the channels. Before continuing, however, let us briefly 
describe a general method to find correctable codes based 
on this tool that may lead to more general applications. 
Let £ = {E a } be an arbitrary quantum channel on B{TL). 
Then all correctable codes for £ may be obtained by ap- 
plying the following procedure: 

(i) For all a, b, find scalars A such that PE\E\,P = XP 
for some projection P. 

(ii) For all A from (i), find the projections P. 

(iii) The set of intersections of ranges of projections P a b 
over all distinct pairs a, b from (ii), corresponds 
precisely to the set of correctable codes for £. 

We emphasize that in QEC, the projections P from (ii) 
must be explicitly identified for practical applications in 
quantum computing. The idea in (iii) will be expanded 
upon in the discussion below. We have stated the prob- 
lem in this form, because it lends itself to consideration as 
a "compression problem" . Of course, this process is still 
somewhat abstract. In particular, step (ii) will typically 
require taking infinitely many intersections. Neverthe- 
less, there may be a way to avoid these infinities. Let us 
describe one simplification of the process. Namely, while 
the family of operators E^Eb, Wa,b, will not be Hermi- 
tian in general, the solving of (i) and (ii) can be reduced 
to the Hermitian case, the case for which we have the 
most general information. 

Proposition 2 The projection P is a solution to Eq. 0) 
for the family {E^Ej, : a, b} if and only if P is a so- 
lution to Eq. 0) for the family of Hermitian operators 
{ElE a ,T+,T- b : a,b}, where 



T 



a b 



ElE h 



4 E a 



(12) 

T- b = i(ElE b -ElE a ). (13) 



Proof. This simply follows from the fact that the operator 
subspace spanned by an operator A and its adjoint A\ 
is also spanned by the Hermitian operators A + A< and 
i(A-A^). ' □ 
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Randomized Unitary Channels. — The class of "random- 
ized unitary channels" [Tij form a large class of physical 
quantum operations. Such a channel has an operator- 
sum representation of the form £(<r) = YlkPkUkcUt, 
where the Uk are unitary and the {pk} determine a prob- 
ability distribution. The error models £ = {Uk} com- 
monly arise in quantum error correction, quantum com- 
munication and cryptography, quantum information pro- 
cessing and theory, etc, and are of great physical rele- 
vance. 

Definition 1 A bi-unitary channel (BUC) is a ran- 
domized unitary channel £ = {V, W} on a Hilbert space 
H. with an operator-sum representation consisting of two 
unitaries; so 

£(a)=pVaV i + (l-p)WaW\ Va £ 13(H), (14) 

for a fixed p with < p < 1 . 

We identify correctable qubit codes for all such chan- 
nels, and we solve the error correction problem for BUCs 
on four-dimensional Hilbert space. The approach we use 
is constructive in nature and allows for an explicit identi- 
fication of the correctable codes. We shall use properties 
of the higher-rank numerical range discussed in the pre- 
vious section. 

There are four non-zero equations to consider in Eq. Q 
here. In terms of the matrix E, we must solve the fol- 
lowing matrix equation for projections P and matrices 
A = (Ay): 

( pP qrPV^WPX fX u P X 12 P\ 
\qrPW^VP {l-p)P ) \X21PX22PJ 1 1 ' 

where we have written q = ^fp and r = y/1 — p. For 
any projection P, the An and A22 equations are trivially 
satisfied with An = p and A22 — 1 — p. Further, the 
A12 equation is satisfied if and only if the A21 equation is 
satisfied. Specifically, 

qrPV^WP = XP if and only if qrPW ] VP = A.P(16) 

Thus, we can reduce the entire problem to solving a single 
(normalized) equation of the form 

PUP = XP, (17) 

for A and P, where U is a unitary on H. 

This can also be seen by simply noting that a code is 
correctable for the channel £ — {V, W} precisely when 
it is correctable for £' = {1, V^W}. In fact, for ease of 
presentation, when it is convenient we shall assume the 
channel £ has the latter form; i.e., there is a unitary U 
on H such that 

£{<t) =pa+{l~p)Uo-U\ y<reB{H). (18) 

Thus, in the case that H = C 4 , we must solve Eq. Ijl7|) 
for an arbitrary 4x4 unitary U. This follows from the 
next result. 



Theorem 3 Let U be a unitary on H = C 4 . Then 
we have the following characterizations of the numerical 
ranges for U : 

(i) Ai([/) = W(U) is the subset of the unit disk in C 
given by the convex hull of the eigenvalues for U . 

(ii) A2(U) is non-empty and given as follows. Let 
z k = e i9k , with 9 k € [0,2vr) for k = 1,2,3,4, be 
the eigenvalues for U , ordered so that < Q\ < 
O2<O 3 <04< 2tt. 

(a) // the spectrum of U is non- degenerate, so 
9 k + Oj Vfc ^ j, then A 2 (U) = {A} where 
X is the intersection point in C of the line £13 
through Z\ and Z3, with the line £24 through Z2 
and Zi . 

(b) IfU has three distinct eigenvalues, sayOj = Ok 
for some pair j ^ k but 9j ^ 9i otherwise, then 
A 2 (U) = {z 3 }. 

(c) // U has two distinct eigenvalues, each of mul- 
tiplicity two, say z andw, then A2(U) consists 
of the line segment L = [z, w] joining z and w. 

(d) If U has two distinct eigenvalues, one X = z 
with multiplicity three and the other with mul- 
tiplicity one, then A2(U) — {z}. 

(e) // U has a single eigenvalue X — z, then U is 
the scalar operator U = zlU and A2(U) = {z}. 

(iii) A^lU) is non-empty if and only if A^(U) = {Ao} 
is a singleton set and Ao is an eigenvalue for U of 
geometric multiplicity at least three; 

dim(ker(Z7 - A 1)) > 3. (19) 

(iv) A^U) is non-empty if and only if U is a scalar 
multiple of the identity operator. 

Proof. The structure of the standard numerical range 
Ai(U) = W(U) follows from well-known matrix analysis 
theory 0], and the structures of Ak(U), k = 3,4, follow 
from Proposition ^ an d Theorem [3 The case of interest 
is that of (ii). 

To see (ii), notice that Theorem El can be applied to 
show that A 2 (U) is contained in the claimed set for each 
of the subcases. (See Fig 2 for a depiction of the four 
non-trivial cases (a), (b), (c), and (d).) For the converse 
inclusion in each case, we offer a constructive proof of the 
required rank-2 projections. 

To verify case (a) , it suffices to show that the intersec- 
tion point A of £13 and £24 belongs to A2(U). First solve 
the following equations, 

J A = azi + bz 3 , . 

\ X = cz 2 + dz 4 ' 1 Uj 
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for nonnegative scalars a, 6, c, d > such that a + b = 1 
and c + d = 1. For instance, a and b = 1 — a may be 
obtained via the equation 

a = cos 2 9 a = X ~ z \ (21) 
z\ - z 3 

for some angle 9 a . Then define an orthonormal pair of 
vectors {\(/>i), \<fa)} by 

( \(j)l) = cos 9 a \ipi) + sin 9 a \ip3 ) , 99 x 
\ |^a) = cos O c \ih) + sin 9 C \^) ' [ ' 

Then we define a rank- 2 projection 

P = |0iX0i| + |M02|- (23) 

Observe that 

(Ufalfa) = cos9 a z 1 (ipi\<f> 1 ) +sm9 a z 3 (ip 3 \<l>i) (24) 
= az\ + bz 3 = A. (25) 

Similarly, we have 1<^ , 2) = A. Further, we also have 

(Ufalfo) =0= {U<h\fa). (26) 

It follows that PU P — XP, and hence A belongs to A^iJJ) 
as claimed. This verifies case (a). 

In case (b), the rank-2 projection P = I^XV'jl + 
IV'fcXV'fel can be seen to satisfy PUP = zjP. For case (d), 
without loss of generality assume z% = z 2 = z 3 = z ^ z±. 
Then P = J2 3 j= i IV'jX^I satisfies PUP = zP. This veri- 
fies both case (b) and (d). 

Finally, for case (c), we have say z\ = z 2 = z and 
Z3 = Z4 = w. Let A belong to the line segment I = [w, z], 
and, as above, solve the following equations for a, 6, c, d; 

J A = az\ + bz 3 = az + bw , 
1 A = CZ2 + dz4 — cz + dw 

Then define {\(j>i), \<fa)} precisely as in Eqs. I)21I22|I . 
WithP= |0iX<£i| + 1^2X^2 1, we have PUP = XP. Thus, 
A 2 (J7) coincides with L = [w, z] in this case, and the re- 
sult follows. □ 

The construction of rank-2 projections used in the pre- 
vious proof is typical in a certain sense. Ostensibly this 
follows because, as proved above, these sets have no topo- 
logical interior. They are either finite discrete sets, or 
line segments. We shall provide full details on this point 
in the next section. (In fact, this result together with 
the discussion of the subsequent section can be easily 
adjusted to v erify the higher-rank numerical range con- 
jecture from [l5j for arbitrary normal operators on 4- 
dimensional Hilbert space.) Of course, a further exten- 
sion of Theorem |3 to the case of arbitrary Ak(U) would 
greatly benefit from either the verification of the conjec- 
ture, or, at the least, more substantive information on 
these sets for arbitrary unitary operators. 

We can establish the following as a consequence of the 
previous proof. 




FIG. 2: Numerical ranges Ai and A2 of a unitary operator 
U 6 U(4) and the structure of the spectrum of U: a) generic 
case; b) double degeneracy; c) both eigenvalues doubly de- 
generated; d) triple degeneracy. 



Theorem 4 Let £ = {V, W} be a BUG on a Hilbert 
space H. with dim7i > 4. Then there are 2- dimensional 
code subspaces CofTL such that C is correctable for £. 

Proof. As in Eq. 1)18(1 . we may assume the channel is of 
the form £ = {11, U}. Thus, the theorem will be proved 
if we can show that A2QJ) is non-empty for an arbitrary 
unitary operator U. 

In the case TC — C 4 , this result follows directly from 
case (ii) of Theorem 01 Indeed, observe that A 2 {U) is 
shown to be non-empty in each of the subcases of (ii). 
The analysis of this case may be adapted to the case of 
arbitrary 7i. To see this, note first that if there is degen- 
eracy in the spectrum of U, say the eigenvalue A = Zj 
has multiplicity at least two, then we can simply choose 
a rank-2 sub-projection of the eigen-projection for Zj to 
show that A 2 (U) is non-empty in this case. 

On the other hand, if there is no degeneracy in the 
spectrum of U , then we may find four distinct eigenval- 
ues for U lying on the unit circle, Zj — e l9j , j — 1, 2, 3, 4, 
ordered so that < 9\ < 9 2 < 9 3 < 9 4 < 2ir. As 
in the previous proof, we can now use the constructive 
"eigenvalue pairing" approach to show that the intersec- 
tion point of the line through z\ and z 3 , with the line 
through z 2 and Z4, belongs to the set A 2 {U). This com- 
pletes the proof. □ 

Let us note that Theorem can also be used in a nega- 
tive fashion. For instance, consider the class of two-qubit 
randomized unitary channels with at least three distinct 
unitary error operators. We can conclude from Theo- 
rem |21 that the set of such channels for which there exist 
correctable codes forms a set of measure zero within the 
set of all such channels. (Contrast this with the case of 
BUCs in Theorem H) Indeed, if £ = {1, U, V} is a two- 
qubit channel with U 7^ V, and neither equal to a scalar 



6 



multiple of 11, then there are three pertinent equations to 
solve of the form given in Eq. I|17|l: namely, for the uni- 
taries U, V, and WV. But we know that in the generic 
non-degenerate case, the set A 2 will be a singleton, and 
that there will almost never be a projection of rank at 
least two that simultaneously solves Eq. (fTTfl for all three 
operators. 

Correctable Codes for the Two-Qubit Case. — We first 
discuss a simple example to illustrate the connection with 
the stabilizer formalism. Let £ = {H4, ZZ} be the error 
model on Ti — C 2 <g> C 2 for a channel that leaves states 
alone with some probability, and applies the Pauli phase 
flip operator Z — |0)(0| — 11)(1| on each qubit with another 
probability. We have chosen this particular example to 
simplify the presentation, but a similar analysis applies 
to any of the two-qubit Pauli-error models. 

The stabilizer formalism applied to the error model 
£ = {J4, ZZ} obtains codes by considering the joint 
eigenspace structure for the error operators {l4,ZZ}, 
and hence just ZZ = Z <g) Z in this case. Through this 
approach one can find the rank-2 projections 

P 1 = |00)(00| + |11)(11| (28) 
P_a = |01)(01| + |10)(10|, (29) 

as correctable codes for £ = {I4, ZZ}. Indeed, Eq. (J2J) is 
satisfied since P ±1 ZZP ±1 = ZZP ±1 = ±P ±1 . 

On the other hand, the higher-rank numerical range 
approach captures these projections, and more. Specifi- 
cally, as ZZ is Hermitian, Theorem ^ can be applied to 
obtain h.i(ZZ~) = [—1,1]. Thus, there is a family of rank- 
2 projections, each of which yields a correctable qubit 
code for £ — {14, ZZ}, for every element A G [—1, 1]. For 
instance, given a, 6, c, d > such that a + b= l = c + d, 
we can define 

= v^|00) + Vb\01) (30) 
= >/3|ll)+Vd|10), (31) 

and define P = \ipi){ipi\ + |^ 2 )(M Then PZZP = XP 
(though note it is not true in general that PZZP = 
ZZP), and it follows that PTt = span{|V>i), is cor- 
rectable for £ — {i4,ZZ}. Notice that Pi is captured 
by b = d = and P_i by a = c = 0. However, the 
cases a b and c ^ d are not obtained as stabilizers of 
two-qubit Pauli-error models. 

Remark 1 There are more projections that yield cor- 
rectable codes for £, and we shall indicate this in the 
context of the next example. We first wish to emphasize a 
basic point exposed by this example. The operator equa- 
tion PAP = XP does not necessarily imply AP = XP; 
i.e., if the compression of an operator A by a projection 
P is a scalar multiple of P, then not necessarily is the 
restriction of A also a scalar multiple of P. In terms 
of the block matrix decomposition of A with respect to 



H = PH ® P these two statements can be phrased 
visually as follows: 



PAP = XP if and only if A 



AP = XP if and only if A = 




It is evident that Eq (|32|l allows for more possibilities. 

To be more precise, given an error model £ = {E a }, 
one way to seek correctable codes for £ is to consider the 
joint eigenspaces for the operators E a (as is done in the 
stabilizer formalism applied for Pauli errors); i.e., pro- 
jections P such that E a P = X a P for some scalar A a , for 
all a. (Recall also that by Proposition [21 the E a may 
be assumed to be Hermitian.) Then necessarily Eq. 
is satisfied for all a, b. However, PE\E^P = X a bP does 
not in general imply that P is a joint eigen-projection 
for E a and Eb. Thus, the previous example, and those 
that follow, are illustrative of a more general phenom- 
ena. Namely, to capture all possible correctable codes, 
we need to consider "compressions" of operators instead 
of restricting ourselves to "restrictions" . 

Let us discuss, in more detail than the previous ex- 
ample, the correctable code structure for the Pauli-error 
model given by £ — {11 4 , Z{\ on H — C 2 <E> C 2 , where 

Z 1 =Z®1 3 = |00)(00| + |01)(01| - |10)(10| - |11)(11|. 

Here the structure of A 2 (Zi) = [—1,1] determines the 
code structure for £. Let us consider the case A = € 
A2(Zi). Let P+ (respectively P_) be the projection onto 
the eigenspace V + = span{|00), 1 01) } (respectively V_ = 
span{|10), |11)}). Then a typical rank-2 code for £ is 
given as follows. Let {|V , ±)i|0±)} be an orthonorma! 
pair of vectors in V± . Define 

C = span{|^+) + |V-), |^+) + \4>-)}. (34) 

One can verify that PZ\P = = OP. For instance, if 
= |^ + ) + |V>->, then 

(Z^IV) H^+l^+)-(^-IV>-) =0. (35) 

On the other hand, suppose C is a two-dimensional 
subspace, with projection P, such that PZ\P = 0. Then 
it follows that C = PH = P+PH® P-PH, and that each 
of the subspaces P±PTL is two-dimensional. Hence, as 
two-dimensional subspaces of V±, we have P±PTi = V±. 
Thus, C necessarily has the form given in Eq. J3"4"fl. (See 
|l5| for a more complete discussion on this point in the 
arbitrary Hermitian case.) 

We finish by characterizing the correctable codes for a 
generic BUC on two-qubit Hilbert space. As above, for 
error correction purposes, we may assume such a channel 
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has the form £ = {I4, U}, where U is a unitary on Ti. = 
C 4 . In the generic case, the spectrum of U will be non- 
degenerate, and so we can assume the eigenvalues for U 
are four distinct complex numbers {z\, Z2, z 3 , 24} , ordered 
so that < argzi < ... < argz 4 < 2tt. Let be 
corresponding eigenvectors; U\ipj) — Zj\ipj). 

By Theorem |3 the correctable qubit codes for £ cor- 
respond to the projections P that satisfy Eq. Ijl7(l for the 
unique A that belongs to A 2 (U). Without loss of gen- 
erality, we shall assume A = 0. (The reader will notice 
that the following argument applies to normal operators. 
Hence, if A 2 (U) = {A}, then we could replace U by the 
normal operator U — AII4.) 

Given that A2(f/) = {0}, we claim that 



p = 16X61 + 16X61 

is a rank-2 projection such that 
PUP = 



(36) 



(37) 



if and only if there are angles a k , /3 k , Ik, #fcj, for fc = 1,2 
and j — 2, 3, 4, such that 



16) = a k \ifo) + b k \^ 2 ) + c k \ip 3 ) + 4|V>4>, (38) 



where 

a k = cos a k cos f3k 

bk = e l8k2 sina fc cos7fc 

Cfe = e l6k3 cosa fc sm/3 k ' 

d k = e l0ki sin a k sin j k 

and the following equations are satisfied; 



cos 2 f3 k zi + sin 2 (3 k z 3 = 
cos 2 7fcZ2 + sin 2 7fcZ 4 = 



(39) 



and 



aia 2 zi + bxb 2 z 2 + ac 2 z 3 + did 2 z 4 = 
a x a 2 zi + b 4 b 2 z 2 + o[c 2 z 3 + d x d 2 z 4 = 



(40) 



(41) 



To verify the sufficiency of these constraints for 
Eq. (PZJ, it must be shown that (Z76I6) = for k, I = 1, 2 
when £1 and 6 are defined as in Eq. (|38|l and the sub- 
sequent equations. This follows from Eqs. (|40I41|) . and 
we leave this computation to the interested reader. Note 
the special case of this construction in which the vectors 
are obtained from the eigenvalue-pairing construction as 
in Eq. (|2*2*|l . (In that case, Eq. l|4*T)> is trivially satisfied.) 

On the other hand, for necessity, suppose P and 6)6 
are given as in Eq. (|3*5|l and Eq. (J2HJ, and that PUP = 
0. Then Eq. (|41|l comes as a direct consequence of the 
identity (£7616) = = (17616) • Further, for k = 1,2 
we have 

= (E/616) (42) 
= \a k \ 2 z 1 + \b k \ 2 z 2 + \c k \ 2 z 3 + \d k \ 2 z 4 . (43) 



In particular, this implies that 

\a k \ 2 z 1 + \c k \ 2 z 3 = -\b k \ 2 z 2 ~ \d k \ 2 Zi = 0, (44) 

as this equation describes the intersection point (A = 0) 
of the line through zi,z 3 with the line through z 2 ,z 4 . 
This yields Eq. (|40|l . and for succinctness we shall 
leave the verification of the specific form of the scalars 
a k ,b k ,c k , d k given by Eq. I|39|l to the reader. 

Let us briefly summarize how one can find all cor- 
rectable codes for any given bi-unitary channel £ = 
{V, W} on H = C 4 : 

(i) Compute the set of compression- values A2(^VT) 
via Theorem |3 

(ii) For each compression- value A from (i), the family 
of projections P that satisfy PV^WP = XP may 
be obtained as in the discussion of this section. 

{Hi) The subspaces corresponding to ranges of projec- 
tions from (ii) axe precisely the correctable codes 
for£ = {V,W}. 

Conclusion. - We have solved the fundamental er- 
ror correction problem in quantum computing for bi- 
unitary channels on two-qubit Hilbert space. This was 
accomplished through an application of a new tool - the 
"higher-rank numerical range" - that has been recently 
developed to solve algebraic compression problems. We 
have shown that, in the case of Pauli-error models, this 
approach captures codes not obtained through the sta- 
bilizer formalism for Pauli matrices. We also discussed 
further applications to more general quantum channels 
on larger Hilbert spaces. As an example, we constructed 
qubit codes for bi-unitary channels on arbitrary Hilbert 
spaces. (Compression error correcting codes for bi- 
unitary channels on arbitrary Hilbert spaces will be an- 
alyzed in a forthcoming article 0.) We also discussed 
how this approach can be used to establish negative re- 
sults. 

To apply the information we have derived for the 
higher-rank numerical ranges in the Hermitian case more 
generally, a better understanding of the intersections of 
projections in Eq. @ is required. As another avenue to 
more general applications, a complete understanding of 
the higher-rank numerical range for the case of normal 
operators, or even unitary operators, could help in spe- 
cial cases such as the class of randomized unitary chan- 
nels. We have not considered possible implications of the 
higher-rank numerical ranges to problems in the new pro- 
tocol for error correction called "operator quantum error 
correction" [jfl El| . In this scheme, a characterization 
of correction has been obtained |2lj, |22| that generalizes 
Eq. J2|, and it may be possible to apply these tools to 
that more general setting. 
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